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Abstract 

The constraints of fixed-target DIS data in fits of parton distributions including QCD corrections 
to next-to-next-to leading order are studied. We point out a potential problem in the analysis of 
the NMC data which can lead to inconsistencies in the extracted value for as{Mz) and the gluon 
distribution at higher orders in QCD. The implications for predictions of rates for Standard Model 
Higgs boson production at hadron colliders are investigated. We conclude that the current range of 
excluded Higgs boson masses at the Tevatron appears to be much too large. 
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The Higgs boson is the last missing cornerstone of the Standard Model (SM). Searches for 
the Higgs boson are in the very center of the experimental activity at the current hadron colliders. 
At the moment, from the combined data of the Tevatron experiments exclusion limits for the SM 
Higgs boson are derived in a certain mass range UJ, while the LHC experiments are in the process 
of improving their discovery (or exclusion) potential with increasing integrated luminosity. At 
the Tevatron and the LHC the Higgs boson can be produced in a large variety of channels, with 
the gluon-gluon fusion process dominating by roughly one order of magnitude over vector-boson 
fusion or Higgs-Strahlung. Precision predictions for the respective production rates are a key 
ingredient in the experimental searches, the higher-order radiative corrections usually being known 
to next-to-next-to-leading order (NNLO) in QCD and to next-to-leading order (NLO) as far as 
electro-weak corrections are concerned (see e.g. [|2l|3l). As favorable features, predictions based on 
higher-order quantum corrections display an apparent convergence of the perturbative expansion 
and a substantially reduced dependence on the choice of the factorization and renormalization 
scales. For the particular case of Higgs boson production in gluon-gluon fusion even the NNLO 
corrections in QCD are still sizable, e.g. roughly 30% for the total cross section, so that NNLO 
accuracy [|4]-|6l is mandatory. 

Phenomenology at hadron colliders, however, also has to address the uncertainty due to the 
non-perturbative parameters, such as the parton distributions (PDFs), the value of the strong cou- 
pling constant UsiMz) and the mass m of the heavy quarks. It has become obvious, that currently 
the largest differences between the various predictions of the Higgs boson cross sections at Teva- 
tron and the LHC are of precisely this origin [|7l[8]|. In this Letter we investigate this point in 
detail. We are concerned here with the value of the strong coupling constant as(Mz) and the PDFs 
as determined in global fits of PDFs and we would like to pin down the source of the resulting 
differences between the PDF sets of ABKM [Q'.fTOl and others. The issue of precision input for 
the value for heavy-quark masses m has recently been solved by using the running mass in the 
MS scheme ifTH . 

PDFs as determined in global fits rely on a variety of data predominantly from deep-inelastic 
scattering (DIS) experiments in order to cover the entire kinematic range in the parton momentum 
fractions Xp. Global PDF fits also combine scattering data with different beams and different 
targets to allow for the separation of the individual quark flavors. Current Higgs boson searches 
probe the PDFs at scales fi of the order of the typical values of the Higgs boson mass Mfj, e.g. say 
// = 165 GeV, and in an effective x^-range determined by (xp) = Mfj/ yfs, where is the center- 
of-mass energy of the collider. The production region at Tevatron is governed by average values 
of (xp) ~ 0.1, while those at e.g. (xp) ~ 0.03 are characteristic for the LHC at = 7 TeV. In this 
Xp-range the relevant experimental constraints on the PDFs are to a great extent due to DIS fixed- 
target experiments (BCDMS UMM, SLAC UM, NMC O, etc). Thus, the processing of these 
data in global fits as well as any assumptions being made must come under scrutiny. As a matter 
of fact, as will be shown below, differences in the treatment of higher-order radiative corrections 
to fixed-target DIS data can be made responsible for the bulk of the deviations in cross section 
predictions based either on the PDF set ABKM MM or on MSTW ^B, the latter being the basis 
of the current Tevatron Higgs searches [[U. 

The fixed-target data are typically provided as differential cross sections of charged-lepton DIS 
off nucleons. In the neutral current case the latter can be written in the one-photon exchange 
approximation as. 
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Figure 1 : The ratio of cross sections R = (Tl/ctt for longitudinally to transversely polarized virtual photon 
DIS as a function of for different values of x. /?nmc denotes the results of lil5J and /?abkm the QCD 
computation to the order indicated. The dashed line corresponds to the result of ||9l. 



where a is the fine structure constant, the (space-like) four-momentum transfer squared, M 
the proton mass and m/ the mass of the incident charged lepton. The Bjorken scaling variable 
is denoted by x and the inelasticity as y (see e.g. [il7i ). The differential cross section in Eq. ([T]) 
depends on the DIS structure functions F2 and Fl. The dependence on Fi can also be parametrized 
by the ratio of the longitudinally to transversely polarized virtual photon absorption cross sections, 
R = (Tl/ctt- The perturbative expansion for the DIS structure functions F2 and F^ in QCD reads 

00 00 
F2 = J]^jFf, F, = a.£a.'Ff, (2) 

/=o 1=0 

with the higher-order corrections being known to NNLO for the PDF evolution [I181I19L as well as 
for the Wilson coefficients of F2 [I20ti22l and Fl [l23l|2l (see also and Refs. therein). Since 
we investigate Higgs boson production to NNLO in QCD, a consistent treatment of the PDFs and 
of the fixed-target DIS data therefore also requires the NNLO corrections for both F2 and F^. 

There exist two possibilities for including fixed-target DIS measurements in global PDF de- 
terminations. One consists of using the differential cross section d^cr/dxdQ^ (of course with all 
electro-weak corrections applied, as required by the respective set of experimental data), i.e. the 
left hand side of Eq. ([T]). This is the procedure of ABKM [9]. Alternatively, one may work directly 
with the published values for F2 extracted from the data for the cross section using the right hand 
side of Eq. ([T]). Although formally equivalent, there are important differences between the two 
approaches, if the latter one does not account on equal terms for the higher-order QCD corrections 

'Note that the correct (9(a'i^) result for was only available after the final publication of the BCDMS data II12I13I . 
but before the NMC analysis flSl . 
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to F2 and F^. This may lead to a significant inconsistency in the PDF fit in particular in the case of 
NMC data analysis. To clarify this point, let us briefly recall a few essentials. NMC was a muon 
beam experiment at CERN with beam energies of 90, 120, 200, and 280 GeV and its data fills 
the gap in the {x, 2^)-kinematics between the SLAC [IT4l| and the HERA measurements at scales 

< 10 GeV^. As such it provides a valuable constraint on the gluon PDF at x > 0.001. 

The extraction of R (or Fi) needs at least two cross section measurements at different beam 
energies for a given x and in order to determine the longitudinal component of the cross section 
from the dependence on y in Eq. ([T]). For NMC, this is not possible in the full kinematic range, 
because the sensitivity to R is substantial only at large y, which implies small x. Thus, only for 
jc < 0.12 NMC has extracted a value for R from its data, while for jc> 0.12 almost all NMC data 
are at j < 0.40 with little sensitivity to R. In this region, x > 0.12, NMC has taken i?i99o from 
lfT4ll . which is based on an empirical parameterization of the SLAC data motivated by QCD and 
including higher-twist terms at large x (see also Il26ll27 ^). A second important issue is concerned 
with the accuracy of QCD perturbation theory. As a matter of fact, the values of R determined by 
NMC rely on leading-order (LO) QCD predictions only. However, since several years two more 
orders in perturbation theory are known for F2 and Fl [|24ll25ll and these higher-order Wilson 
coefficients contain the non-trivial dependence and induce big corrections in particular at small 

In Fig. \T\ we display results for R as published by NMC in comparison to values calculated 
with the ABKM PDFs ^ in various orders of perturbation theory. The two x values, x = 0.11 and 
X = 0.14, are chosen to illustrate the different analysis strategies of NMC, i.e. either an extraction 
of R from its data (x = 0.1 1) or the use of i?i99o from yjl (x = 0.14). We find that the value of R 
obtained with the ABKM [|91 fit is in good agreement withi?i99o atx> 0.12. This is because similar 
sets of data (SLAC f\M) are used in both fits. Fig. [T] also illustrates the impact of the higher-order 
Wilson coefficients, i.e. F^^^ and F^^^ from Eq. Q, which lead to good perturbative stability of 
the NNLO prediction even at small values of Q^. However, the value of R^mc at jc < 0.12 is quite 
diff"erent from Riggo at jc> 0. 12. In contrast to the R values computed with ABKM [9], ^nmc does 
not depend on at x < 0.12, see Fig. [T] (left), due to assumptions made in the NMC analysis. 
Thus, it is obvious that high precision PDF fits to NLO or NNLO in QCD need to be based on the 
NMC data for the differential cross sections d^cr/dxclQ^ rather than on the NMC results for F2 and 
R, because the latter approach is simply inconsistent. Of course, a similar statement also holds for 
the analysis of other fixed-target data, where modern parametrization of R with higher-order QCD 
corrections have to be applied as well and, e.g. in the case of BCDMS [|T2l[T3l , have an impact on 
the valence quark PDFs [i28il . 

It is interesting to investigate the consequences of the two alternative treatments of the NMC 
data, which covers the range of jc ~ 0.001 ...0.1, and thus is of great importance for the Higgs 
boson production at current hadron colliders. For a quantitative analysis we perform a variant 
of the ABKM fit BUl with the NMC data for the cross section replaced by the data for F2. The 
results for the values of the strong coupling constant as{Mz) are presented in Tab. [IJ Interestingly, 
the (inconsistent) NNLO variant of the ABKM fit based on the NMC data on F2 and R yields a 
value of Q'i(Mz) = 0.1170, bigger than the default ABKM value by -1-0.0035 and rather close to 
to MSTW [16| and the present world average of OsiMz) [29| El- The resulting shift corresponds 
to more than +2cr standard deviations. In the NLO case, the shift is smaller, -1-0.0009, whereas it 

^Note that values of a^iMz) from NLO, NNLO and N^LO determinations contribute to this average. In world 
analyses of DIS and other hard scattering data the extracted NNLO values for as{Mz) are systematically lower than 
the con-esponding NLO ones, e.g. by -0.0044 (El, -0.0031 Hg) or -0.0014 1281. 
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ff.v(Mz) 


a,.(Mz) with ctnmc 


a,{Mz) with F^^c 


difference 


NLO 


0.1179(16) 


0.1195(17) 


+0.0026 ^ 1.5cr 


NNLO 


0.1135(14) 


0.1170(15) 


+0.0035 =^ 13(T 


NNLO +Fl at 


0.1122(14) 


0.1171(14) 


+0.0050 =^ 3.6cr 



Table 1 : The values of the strong coupling as{Mz) obtained in global fits of PDFs from variants of the 
ABKM fit [9]. The order of perturbation theory is indicated in the left column and in the two central ones 
the treatment of the NMC data 1 15 ] , i.e. a fit to the measured cross sections or to the DIS structure function. 
The right column gives the absolute difference and the relative one in terms of standard deviations. The 
value in bold corresponds to the published result in 191 . 




\i=\6S GeV 



Figure 2: The absolute uncertainty of the gluon PDF as a function of x at the scales // = 2 GeV (left) and 
li-\65 GeV (right) for the ABKM fit of Ref. 121 (shaded area) compared to the difference with the variant 
of ABKM fit with used (dashes). 



becomes even larger if the 0{as^) corrections for the Wilson coefficients of Fi are included. Then, 
the difference between the two treatments amounts to +3.6cr. This is to be expected, because 
the alternative treatments of the NMC data are almost equivalent at LO (see Fig. [U and deviate 
more and more as we include higher and higher orders for Ft. The as(Mz) values resulting from a 
consistent treatment of the NMC data (left column in Tab.[T]) are in full agreement with other recent 
high precision determinations 1301 [STl . It should also be mentioned here, that the values ofx^ in 
all variants of the fit are very similar. The variations are roughly +10 units, which is statistically 
insignificant given the large number of data points in the fit. This means that the variation in the 
ansatz is fully compensated by the changes in the PDFs and the value of OsiMz). 

In Fig. [21 we plot the change in the gluon PDF G(x) due to the choice of the NMC data represen- 
tation. The variant with a fit to F^^'^ displays significant deviations. At the initial scale // = 2 GeV 
(Fig. [21 left) it effectively leads to a larger gluon in the range .x: > 0.1. One should keep in mind 
here, that the gluon PDF at larger scales is actually sensitive to all values of x larger than 0.1, 
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^=165 GeV 




Figure 3: The gluon luminosity Lgg = gi^g (weighted by a factor a/) at the scale /j. = 165 GeV as a 
function of the Higgs boson's rapidity Y. The solid line denotes the result of Ref. ||9l and the dashed line the 
variant with a the fit to F^'^*^. 



because the physical observables emerge as convolutions with the respective Wilson coefficients. 
Due to the QCD evolution from /i = 2 GeV to the scale /i = 165 GeV this excess of the gluon PDF 
then extends to even smaller x > 0.05 (Fig. [2l right), which is precisely the range in x relevant for 
Higgs production at Tevatron and the LHC. Thus, in the inconsistent variant of the fit to NMC 
data, one obtains both, a larger value of as{Mz) and a larger gluon PDF. This matches with the 
observed differences between the PDF sets currently available at NNLO in QCD. For the gluon 
PDF in the relevant x range, .jc ^ 0.1, (see e.g. Fig. 2 in Ref. [8J), we find good agreement between 
ABKM and the results of HERAPDF [32]. The latter are obtained from a fit without NMC data. 
On the other hand, no agreement exists with JR [[331 and MSTW [[T6l . These fits both use the 
NMC results for F2. Remarkably, in comparison to ABKM, both the gluon PDF and the a? value 
of MSTW are larger. In this context it should also be stressed the initial conditions for the gluon 
PDF are significantly correlated with the value of UsiMz) which determines the speed of the QCD 
evolution. Especially at large x a strong anti-correlation is observed (see e.g. Tab. 2 in Ref. [(9|[), 
so that a smaller Qr^(Mz) value implies a larger gluon PDF and vice versa. 

Finally, we would like to summarize the impact of the different variants to treat the NMC 
data on the predicted Higgs boson cross sections at the Tevatron and the LHC. We focus on the 
dominant channel through gluon-gluon fusion and illustrate the cumulative effect of a larger gluon 
PDF and a larger Us value. The Bom contribution in gluon-gluon fusion is proportional to and 
the gluon luminosity Lgg = G®G. In Fig. [3] we plot the product cus^Lgg at the scale = 165 GeV 
as a function of the Higgs boson's rapidity Y. The difference between the ABKM prediction and 
the (inconsistent) variant with a fit to F^MC ^^^^ amounts to an increase of roughly 20% at central 
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a{H) 


cr{H) with ctnmc 


ct{H) with f NMC 


difference 


NLO 


0.206(17) pb 


0.225(18) pb 


0.019 pb ^ l.lcr 


NNLO 


0.253(22) pb 


0.309(24) pb 


0.056 pb ^ 13o- 


NNLO +Fl at 0(a,^) 


0.242(22) pb 


0.310(24) pb 


0.068 pb - 2.8cr 



Table 2: The predicted cross sections for Higgs boson production in gluon-gluon fusion with Mu - 
165 GeV at Tevatron ( = 1.96 TeV) obtained with the PDFs from variants of the ABKM fit |9]. The 
order of perturbation theory is indicated in the left column and in the two central ones the treatment of the 
NMC data |[T5l . i.e. a fit to the measured cross sections or to the DIS structure function. The right col- 
umn gives the absolute difference and the relative one in terms of standard deviations. The value in bold 
corresponds to the published result |[8||. 



CT{H) 


(t{H) with (Tnmc 


cr{H) with F^mc 


difference 


NLO 


5.73(17) pb 


5.95(18) pb 


0.18 pb ^ l.Oo- 


NNLO 


7.05(23) pb 


7.70(23) pb 


0.65 pb ^ l.lcr 


NNLO +Fl at 0(Qf,^) 


6.84(21) pb 


7.68(23) pb 


0.84 pb ^ 3.7o- 



Table 3: Same as Tab.[2]for the LHC ( - 7 TeV). 

rapidities. In order to quantify this enhancement for total cross section predictions we present in 
Tabs. [21 and [3] the respective numbers. For the Tevatron (Tab.O the NNLO QCD prediction based 
on the fit to the F^MC ^^^^ yields a cross section value which is 22% larger than the one from 
ABKM C9|. This corresponds to a shift of +2.3cr standard deviations of the combined uncertainty 
on Us and the PDFs, a difference which still increases slightly if the Wilson coefficients for Fl at 
0{as^) are included. At NLO however, the difference is of the order of +lcr only, which is in line 
with the previous observations in the determination of (Mz), cf . Tab. [Tl For the LHC (Tab. [3]) at 
^fs = l TeV center-of-mass energy the same pattern emerges. At NNLO the inconsistent treatment 
of the NMC data in the fit leads to a cross section which is 9% larger than the ABKM prediction H 
and the difference amounts to -l-2.7cr standard deviations. 

The results in Tabs. |2] and [3] provide a potential explanation for the significant spread in the 
predicted Higgs cross sections, especially between the ABKM and MSTW PDF sets, where the 
differences are largest. In the present exclusion region for Higgs masses around Mh = 165 GeV, 
the MSTW prediction at NNLO in QCD for Tevatron is +35% higher than the one of ABKM, i.e. 
a +4. Oct deviation in the combined as and PDF uncertainty. At the LHC with ^fs = 1 TeV the 
respective MSTW prediction is still +12% higher, which corresponds to a deviation of +3.6(T (see 
the detailed study in [8] for numbers). The different handling of the NMC data in the global fits 
that is characterized by either accounting for or neglecting higher-order corrections to Fi accounts 
for the bulk of the observed deviations. 

In summary, we have highlighted the importance of fixed-target DIS data for predictions of 
rates for SM Higgs boson production at hadron colliders. The use of the NMC data in global 
PDF analyses allows for different choices. It is preferable to rely on the differential cross sections 

as we have shown that a direct fit to F2 from NMC leads to inconsistencies 
at higher orders in QCD. We have illustrated the implications of these options for the determination 
of as(Mz) and PDFs in a global fits and we have computed the rates for Higgs boson production 
at Tevatron and LHC with the results of these fits. The observed differences in the Higgs cross 
section allow us to understand the deviations in the predictions between the currently available 
NNLO PDF sets, most prominently between ABKM and MSTW, which differ by roughly 4cr in 
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the combined uncertainty for the PDF parameters and a 5. The details of the DIS fixed-target data 
analysis are therefore an important issue for the interpretation of the Tevatron data and for limits 
on the mass of a SM Higgs boson [IJ. The current range of excluded Higgs boson masses appears 
to be much too large. It could easily be overestimated by a factor of two based on the reduced rate 
for the Higgs boson signal alone and this topic needs urgently further investigation. Potentially, 
studies of the projected sensitivities for SM Higgs production at the LHC are also afl'ected by this 
concern. In any case, it will be mandatory to base upcoming SM Higgs searches at the LHC on 
parton luminosities from consistent global PDF fits. 
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